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I. INTRODUCTION 

Interest in the temporal and spatial variation of major 
constants of physics has been recently revived by astro- 
nomical data which seem to suggest a variation of the 
electromagnetic constant a = /he at the 10~^ level for 
the time scale 10 billion years, see Q (a discussion of 
other limits can be found in the review and references 
therein). However, an independent experimental confir- 
mation is needed. 

The hypothetical unification of all interactions implies 
that variation of the electromagnetic interaction constant 
a should be accompanied by the variation of masses and 
the strong interaction constant. Specific predictions need 
a model. For example, the grand unification model dis- 
cussed in Q predicts that the quantum chromodynamic 
(QCD) scale Aqcd (defined as the position of the Landau 
pole in the logarithm for the running strong coupling con- 
stant) is modified as follows: SAqcd/^qcd ~ SASa/a. 
The variation of quark and electron masses in this model 
is given by Sm/m ~ 706a/a. This gives an estimate for 
the variation of the dimcnsionless ratio 



{m/AqcD) a 



(1) 



This result is strongly model-dependent. However, the 
large coefficients in these expressions are generic for 
grand unification models, in which modifications come 
from high energy scales: they appear because the running 
strong coupling constant and Higgs constants (related to 
mass) run faster than a. This means that if these models 
are correct the variation of masses and strong interaction 
may be easier to detect than the variation of a. 

One can measure only variation of the dimcnsionless 
quantities, therefore we want to extract from the mea- 
surements variation of the dimensionless ratio mq/ Aqcd 
where ruq is the quark mass (with the dependence on the 
normalization point removed). A number of limits on 
variation of mq/ Aqcd have been obtained recently from 
consideration of Big Bang Nucleosynthesis, quasar ab- 
sorption spectra and Oklo natural nuclear reactor which 
was active about 1.8 billion years ago 0, IE B (see 
also 

II IS E in II3)- Below we consider the limits 



which follow from laboratory atomic clock comparison. 
Laboratory limits with a time base about a year are es- 
pecially sensitive to oscillatory variation of fundamental 
constants. A number of relevant measurements have been 
performed already and even larger number have been 
started or planned. The increase in precision is very fast. 

It has been pointed out by Karshcnboim 13] that mea- 
surements of ratio of hyperfine structure intervals in dif- 
ferent atoms are sensitive to variation of nuclear mag- 
netic moments. First rough estimates of the dependence 
of nuclear magnetic moments on mq/ Aqcd and limits 
on time variation of this ratio have been obtained in our 
paper Using H, Cs and Hg+ measurements |T^ ITU - 
we obtained the limit on variation of mq/ Aqcd about 
5 • 10"^'^ per year. Below we calculate the dependence 
of nuclear magnetic moments on mq/ Aqcd and obtain 
the limits from recent atomic clock experiments with hy- 
perfine transitions in H, Rb, Cs,Yb+,Hg+ and optical 
transition in Hg+. It is convenient to assume that the 
strong interaction scale Aqcd does not vary, so we will 
speak about variation of masses. We shall restore Aqcd 
in final answers. 

The hyperfine structure constant can be presented in 
the following form 



const X 



[a^ Frei{Za)][^—] 



(2) 



The factor in the first bracket is an atomic unit of en- 
ergy. The second "electromagnetic" bracket determines 
the dependence on a. An approximate expression for the 
relativistic correction factor (Casimir factor) for s-wave 
electron is the following 



7(47' - 1) 



(3) 



where 7 — \J\ — (Za)^, Z is the nuclear charge. Varia- 
tion of a leads to the following variation of Frei 0| : 



SFr, 



K 



5a 



K = 



(Za)2(127^ - 1) 



(4) 



(5) 



72(472 - 1) 

More accurate numerical many-body calculations of 
the dependence of the hyperfine structure on a have 



2 



shown that the coefficient K is sUghtly larger than that 
given by this formula. For Cs {Z—55) K— 0.83 (in- 
stead of 0.74), for Rb K^QM (instead of 0.29), for Hg+ 
K=2.2% (instead of 2.18). 

The last bracket in eq. ((SJ contains the dimensionless 
nuclear magnetic moment /i in nuclear magnetons ( the 
nuclear magnetic moment M = u „^^ ), electron mass 

rrie and proton mass mp. We may also include a small 
correction due to the finite nuclear size. However, its 
contribution is insignificant. 

Recent experiments measured time dependence of the 
ratios of hyperfine structure intervals of ^^^Hg"*" and H 
p^, ^^^Cs and ^"^Rb 0| and ratio of optical frequency 
in Hg+ and ^■^^Cs hyperfine frequency |T3|. In the ratio 
of two hyperfine structure constants for different atoms 
time dependence may appear from the ratio of the fac- 
tors Frei (depending on a) and ratio of nuclear magnetic 
moments (depending on mq/ Aqcd)- Magnetic moments 
in a single-particle approximation (one unpaired nucleon) 
are: 



for j = / + l/2. 



M = [gs + (2j - l)gi)/2 



i-g, + (2j + 3)50 



2(J + 1) 



(6) 



(7) 



for j = Z — 1/2. Here the orbital g-factors are gi = 1 
for valence proton and gi = for valence neutron. The 
present values of spin g-factors gs are gp = 5.586 for 
proton and gn = —3.826 for neutron. They depend on 
rUql Kqcd- The light quark masses are only about 1% 
of the nucleon mass {mq = {rnu -f md)/2 « 5 MeV). 
The nucleon magnetic moment remains finite in the chiral 
limit of jUu — md — 0. Therefore, one may think that the 
corrections to g^ due to the finite quark masses are very 
small. However, there is a mechanism which enhances 
quark mass contribution: 7r-meson loop corrections to 
the nucleon magnetic moments which are proportional 
to TT-meson mass ~ yJmqt^QCD\ m7r=140 MeV is not 
so small. 

According to calculation in Ref. dependence of the 
nucleon g-factors on 7r-mcson mass ttItt can be approxi- 
mated by the following equation 



.9(0) 



1 -|- flTOTT -f- bm'i 



(8) 



where a= 1.37/GeV , 6= 0.452/GeV2 for proton and a= 
1.85/GeV , b= 0.271/GeV2 for neutron. This formula 
is a fit of the numerical calculations performed using the 
Cloudy Bag Model (pion field coupled to nucleon quarks). 
This model reproduces leading terms of chiral perturba- 
tion theory for nucleon magnetic moments H^l (small 
limit). The results also agree with lattice calculations 
[lil (performed at - 500 MeV). Eq. © leads to the 
following estimates: 



Sgp 
9p 



-0.174- 



TO, 



-0.087^ 



(9) 



— = -0.213 

gn 



-0.107^ 



(10) 



Eqs. H6I7I9I10() give variation of nuclear magnetic mo- 
ments as functions of mq/ Aqcd- 

The measured time variations of ratios of hyperfine fre- 
quencies depend on two parameters: the ratio of proton 
spin magnetic moment and proton orbital magnetic mo- 
ment (this ratio is proportional to gp) and the ratio of 
proton and neutron spin magnetic moments {Mp/Mn = 
gp/gn)- According to eqs. 110|l the ratio gp/gn practi- 
cally does not depend on mq and seems to be not sensi- 
tive to variation of quark masses and strong interaction. 
However, this conclusion may be misleading. Magnetic 
moments depend also on the strange quark mass m^. In 
the minimal order of the chiral perturbation theory this 
dependence is dominated by mixing p — K. Similar 
process involving meson loop and A does not con- 
tribute to the neutron magnetic moment since K'^ does 
not carry electric charge. As a result lowest order cor- 
rections for proton and neutron magnetic moments are 
very different: for proton a(K)la{T:) = 0.4, for neutron 
a{K)/a{T:) = -0.03 [H 113 (see eq. ® for definition 
of the coefficient a). The proton (or neutron) g- factor 
is a ratio of the proton (or neutron) spin magnetic mo- 
ment to the nuclear magneton efi/2mpC (quantum for 
proton orbital magnetic moment). Therefore, depen- 
dence of gp ~ Mpmp and g„ ~ M^mp on the strange 
quark mass may also appear from the relatively large 
(~ 20%) contribution of TOs to the proton mass mp 
In the case of gp this contribution is of opposite sign to 
that of the lowest order in the chiral perturbation theory 
{a{K)mK) and has comparable magnitude. This cancel- 
lation means that we can not reliably determine depen- 
dence of gp on the strange quark mass m^. In this case 
it would be safer to neglect contribution of the strange 
quark mass to gp. However, we want to keep dependence 
on the strange quark mass in the ratio of neutron and pro- 
ton magnetic moments since in a number of cases this is 
the only effect which provides dependence on fundamen- 
tal masses and strong interaction. Using the lowest or- 
der chiral perturbation theory results [20l l2ll| presented 
above we obtain the dependence on the strange quark 



Q.r- 



-, ~ u.-L . This result should possibly be 

treated as an order of magnitude estimate since K- me- 
son mass is not as small as 7r-meson mass and accuracy 
of the chiral theory may be low in this case. Thus, we 
arrive to the following equations: 



^ ^ -0.087^ 



9v 



ma 



Sg. 



SjTt 

— - -0.107( ^ 

gn mq 



5m, 



-) 



(11) 



(12) 



Note that transferring the entire TO^-dependence into the 
neutron g-factor is mostly a matter of convenience, all 
that we can say is that gn/ gp depends on m^/ Aqcd- Now 
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we can find variation of nuclear magnetic moments using 
eqs. H6I7|I . For all even Z nuclei with valence neutron 
(i99jjg^i7iYb,iiiCd, etc) we obtain ^ = • For i^^Cs 
we have valence proton with j=7/2, I— A and 



^ = 0.22^ 



0.11- 



Snia 



(13) 



For ^'Rb we have valence proton with j=3/2, Z=l and 



5n SrriTr 
— = -0.128 



-0.064^ 
ma 



(14) 



As an intermediate result it is convenient to present 
dependence of the ratio of hyperfine constant A to atomic 

4 

unit of energy E = "^^i (or energy of ls-2s transition in 
hydrogen equal to 3 /8 E) on variation of the fundamental 
constants. We introduce a parameter V defined by the 
realation 



SV_ ^ 5{A/E) 
'V " A/E 



(15) 



We start from the hyperfine structure of "^^^Cs which is 
used as a frequency standard. Using eas. (|2ll3|l we obtain 

VC^^Cs)^a^-'\^r'\^) (16) 

i^QCD J^^QCD 

Here we have taken into account that the proton mass 
TOp oc Aqcd- For hyperfine transition frequencies in 
other atoms we obtain 

VrRb) = a'-'\^r°-°'{^) (17) 



QCD 



'A 



QCD 



VCH)^a\^)-^-^\-^) (18) 



VCH)^a\-^)-'-'\-^)-'-^\^) (19) 
J^QCD Aqcd Aqcd 



VC'''>Hg+) = a^-^{—f-^\^^) (20) 
mq Aqcd 



VC'^Yb+)=a^-'{'^f-^\-^) (21) 
rriq Aqcd 



VC''Cd+) = a^-<^{^f-'\-^) (22) 
niq Aqcd 

Now we can use these results to find limits on variation 
of the fundamental constants from the measurements of 
the time dependence of hyperfine structure intervals. The 
dependence of ratio of frequencies A{^^^Cs)/A{^'^Rh) can 
be presented in the following form 

X{CsRb) = J3St = «"-''K/Aqcd]°-" (23) 



V{Rb) 



Therefore, the result of the measurement [13 may be 
presented as a limit on variation of the parameter X: 



1 dX{CsRb) 
X(CsRb) Jt 



(0.2±7) X lO^'^/year (24) 



Note that if the relation is correct, variation 
of X{CsRb) would be dominated by variation of 
ql Aqcd]- The relation iQJ would give X{CsRb) oc a^. 



For A(i33Cs)/A(H) we have 
X{CsH) 



^^^'^ -a'-''[mq/AQCDr (25) 



ViH) 

Therefore, the result of the measurements [Tsl may be 
presented as 



1 dX{CsH) 



'X{CsH) dt 
For A(i99Hg)/yl(H) we have 



< 5.5 X W'^'^/year (26) 



X{HgH) 



^^"'^ -a^-^K/AQCD]"-^ (27) 



V{H) 

The result of measurement |l4j may be presented as 
1 dX{HgH) 



'X(HgH) dt 



-\<8x 10-^^ /year (28) 



In Ref. [13 the limit on variation of the ratio of hy- 
perfine transition frequencies ^^^Yb'^/^'^^Cs has been ob- 
tained (this limit is based on measurements [l^). Us- 
ing eqs. (|16I21|) we can present the result as a limit on 
X{YbCs) = «0-7K/Aqcd]"-": 



1 dX{YbCs) 
X{YbCs) Jt 



1(2) X IQ-'-^/year (29) 



The optical clock transition energy E{Hg) (A=282 nm) 
in Hg"*" ion can be presented in the following form: 

4 

E{Hg) = const x [-A^]Frei{Za) (30) 

Numerical calculation of the relative variation of E{Hg) 
has given p^ : 



SEjHg) _ Sa 
E{Hg) ■ a 



— ^,-3.2 



(31) 



Variation of the 



This corresponds to V{HgOpt) 
ratio of the Cs hyperfine splitting A{Cs) to this optical 
transition energy is described by 

^^Opt) = ^/j^'^^, = a6K/AQCD]°-"K/AQCD] 



V{HgOpt) 



(32) 



The work [T^ gives the limit on variation of this param- 
eter: 



1 dXjOp t) , 
'x(Opt) dt 



I < 7 X lO^^Vyear (33) 
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Molecular rotational transitions frequencies are pro- 
portional to mf./mp, therefore one should assume V = 
mel ^QCD- For vibrational molecular transtions V — 
{me/ h^QcoY^"^ ■ Comparison of Cs hyperfine stan- 
dard with SFq molecular vibration frequencies was dis- 
cussed in Ref. ■ In this case X{CsVibrations) = 

a2.8[™JAQCD]0-5[m,/AQcz3]"-^ 

The measurements of hyperfine constant ratio in dif- 
ferent isotopes of the same atom depends on the ratio 
of magnetic moments and is sensitive to ruq/ Aqcd- For 
example, it would be interesting to measure the rate of 
change for hydrogen/deuterium ratio where X{HD) = 

K/Aqcd]-°-°^K/Aqcd]°-'^. 

R. Walsworth suggested to measure the ratio of the 
Zeeman transition frequencies in noble gases which gives 
us time dependence of ratio of nuclear magnetic mo- 
ments. Consider, for example ^^^Xe/'^He. For '^He the 
magnetic moment is very close to that of neutron. For 
other noble gases nuclear magnetic moment is also given 
by valence neutron, however, there are significant many- 
body corrections. For ^^^Xe the valence neutron is in 
state which corresponds to the single-particle value 
of nuclear magnetic moment fj, — fin — —1.913. The 
measured value is = —0.778. The magnetic moment of 
the nucleus changes most efRciently due to the spin-spin 



interaction because valence neutron transfers a part of its 
spin < Sz > to core protons and proton magnetic moment 
is large and has opposite sign. In this approximation fi = 
(1 — -I- bfip. This gives b=0.24 and the ratio of mag- 
netic mom ents Y = ^(i29Xe)//i(3He)« 0.76 -f 0.24gp/5„. 
Using eas. (|llll2|) we obtain an estimate for the varia- 
tion of this ratio SY/Y « OASnis/ms. Therefore, in this 
case one can measure variation of /i(^^^Xe)//i('^He) cor- 
responding to variation of A" = [ms/ Aqcd]'^'^ ■ 

Note that an accuracy of the results presented in this 
paper depends strongly on fundamental constant we are 
studying. The accuracy for the dependence on a is few 
percent. The accuracy for rUq/AQCD is about 30%, it 
is limited by the accuracy of a single-particle approxi- 
mation for nuclear magnetic moments. The results for 
uIs/Aqcd should possibly be treated as order of mag- 
nitude estimates. However, the accuracy here may be 
improved using proper QCD calculation 0|. Relation 
between variation of a and tti/Aqcd has been used 
as an illustration only. 
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